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ABSTRACT 

A system of particles hopping on a line, singly or as merged pairs, and annihilat- 
ing in groups of three on encounters, is solved exactly for certain symmetrical initial 
conditions. The functional form of the density is nearly identical to that found in two- 
body annihilation, and both systems show non-mean- field, ~ instead of ~ 1/t, 
decrease of particle density for large times. 
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Diffusion limited reactions in low dimensions have attracted much theoretical in- 
terest [1-11]. Indeed, fluctuation effects are most profound in low dimensions, implying 
deviations from the mean-field rate-equation behavior. Numerical, analytical, and in 
few cases exact results have been obtained for two-body reactions A + A ^ products, 
and A + B ^ products, in d = 1. Particle input and production (fragmentation) have 
also been allowed for in recent studies [4,8,12]. 

Exact solutions are available only for the simplest reactions. Specifically, for 
the two-body annihilation, A + A ^ inert, the most detailed exact results have been 
obtained by methods also used for the Glauber-dynamics Ising models. The particles A 
are then interfaces on a lattice dual to the Ising spin lattice [5,9,13]; the time evolution 
is continuous. Exact results for cellular-automata type, simultaneous-update, discrete 
time dynamics, have also been derived [11]. For reactions A + A A, possibly 
with particle input and fragmentation, continuum-limit results have been obtained by 
various random- walk arguments [7,12], some of which are not limited [7] to d = 1. 

In this work, exact results will be obtained for a certain three-body reaction 
scheme to be defined shortly. The method of solution involves consideration of 
Zfc-symmetric, fc-state "spins" on a linear lattice. Multistate, Potts-spin d — 1 lattice 
has been considered [14] as a model of mixed pairwise annihilation and coagulation of 
interfaces between the fully ordered (T = 0) Potts phases; the limit k = oo was solved 
by random-walk methods [14]. The approach adopted here is quite different: discrete- 
time and space cellular-automaton dynamics will be defined, and the case k = 3 will 
be shown to have certain simplified properties allowing an exact solution. 

It is convenient to use directly the roots-of- unity complex representation of Z^, 
the numbers exp {27rin/k), with n = 0,1,2, ... ,k — l. Our "spin" variables, (Jj{t), take 
on these k distinct values. Here the time t increases in integer steps, t — 0, 1, 2, . . .. 
The spin variables are at even lattice sites j = 0, ±2, ±4, . . ., for even times, and at 
odd lattice sites for odd times. Let 
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CTj (t) = exp [27rinj (t) / k] 



(1) 



We identify the reacting-diffusing particles as follows. On a dual lattice, i.e., in the 
interstice between two neighboring spins at j — 1 and j + 1, we calculate the number 
aj+i/aj-i = aj+iaj_i, where star denotes complex conjugate. The result is put in 
the form exp {27riNj/k), where A^^ = (n^+i — n^-i, mod A;), is a number in the range 
0, . . . , — 1. We then identify the number of particles at j (at time t) as Nj{t) thus 
calculated. Note that the particles occupy the dual lattice sites with even j for odd 
times, and with odd j for even times t. 

Next, we define the reaction-diffusion dynamics. At each time step t t + 1, 
all the Nj{t) particles at the dual-lattice site j hop as a group to the right, to site 
j -I- 1 of the dual lattice, or to the left, to site j — 1, with equal probability. Thus 
the resulting number of particles hopping to site m can be in the range 0, ... ,2k — 
2. More precisely, this number can be 0, Nm-i, Nm+i, or the sum Nm-i + Nm+i- 
However, we now add the reaction process of instantaneous annihilation of k particles, 
at all sites at which there are more than k particles present after hopping. Thus, in 
those sites m which received Nm-i{t) + Nm+i{t) particles, we define N^it + 1) as 



These dynamics can be put in a convenient linear form in terms of the 
Zfe-symmetric spins. 



Here we assigned independent and uncorrelated random variables Cm(0 to each spin 
at each time step t > 0. Stochastic variables Cm{t) take on values or 1 with equal 
probability. For k = 2 this and related linear dynamical rules lead to a factorizable 
hierarchy of relations for equal-time, p-spin correlation functions [11,15]. The factor- 



{Nm-i{t)+N^+i{t), mod/c). 



+ 1) = [1 - Q{t + 1)] ^j-iit) + Cj{t + 



(2) 



3 



ization property will also apply for /c > 2. 

Essentially, the p-spin correlation functions at fixed time t can be calculated from 
relations which only involve p, {p — 1), {p — 2), . . -spin correlations, but do not couple 
to higher-order correlations. The resulting expression are, however, difficult to treat 
exactly, and all the available results for = 2 have been for 2-spin correlations with 
symmetric (translationally-invariant) initial conditions. In terms of the spin variables, 
the simplest tractable initial conditions are fully random ones: spins (Tj{0) are uncor- 
related at different sites j, and distributed identically at eachj', with (cr-, (0)) = s. Here 
s is a complex number which, for k = 3, must be within (or on the boundary of) the 
triangle defined by 1, exp(27ri/3) and exp(47ri/3) in the complex plane. 

The average (...) for times t > will denote both averaging over the stochastic 
evolution, i.e., over the choices of the variables Cm{t)i and over the distribution of the 
initial values Cj(0). However, the s uncorrelated distribution of the initial spin 
values implies a subtle correlation in the initial distribution of the reacting-diffusing 
particles on the dual lattice, as emphasized for A; = 2 in the Glauber-type approach 
[9]. 

Thus, even though our solution can be easily extended to general s, we will put 
s = 0. In fact, we assume that the reacting-diffusing particles are initially distributed 
as follows for k — 3: at each dual-lattice site, one randomly places 0, 1, or 2 particles 
at time t = 0. The average initial density is p(0) = 1, and such distribution ensures 
that each spin has equal probability to have the three allowed complex values (for 
k — 3; extension to A; 7^ 3 is trivial). 

For times t > the particle density will decrease from p(0) = 1. In fact, the 
mean- field rate equation prediction for /c = 3 can be easily derived, and it is p{t) ^ 1/t 
for large times. We now put k = 3 and derive the exact expression for p{t). The result 
yields p{t) ~ l/y^ for large times, similar to the fluctuation-induced modification of 
the asymptotic density in the case of the two-body annihilation A + A ^ inert, which 
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corresponds to another solvable case, k — 2. 

The case = 3 is special due to the property that the density of particles on the 
dual lattice can be expressed in terms of the 2-spin correlations of the Zk spins for 
k = 3 (and if fact also for k = 2) but not for k > 3. Indeed, if = exp {27ri£/k) is one 
of the k roots of unity constituting the representation of Zk>2, then the appropriate 
phase number, £ = 0, 1, . . . , /c — 1, can be calculated as a linear form a a (2) 9 + (1(3)9* 
only for k = 3. Indeed, three complex coefficients a(i),(2),(3) are just the right number 
to be fixed to get the three values £ = 0,1,2 properly. For k > 3, the expression 
required would involve higher integer powers of 9 and 9* and it would no longer be 
linear. For k = 3 the explicit result is 



e=l + ^Re 



{-V3 + i)g 



(3) 



Recalling our rule of calculating the particle number Nj (t) , we conclude that the 
average particle density at site j (of the dual lattice) at time t is given by 



Pit) = 1 + -^Re [(-^ + i)(a,+i(t)cT;_i(t))] (k = 3) (4) 

where we omitted the index j for p{t) due to translational invariance. Relation (4) ap- 
plies only for = 3 and expresses the particle density in terms of the 2-spin correlation 
function of the spins. 

The dual-lattice sites can be empty, singly-occupied, or pair-occupied. Let Psingie(^) 
denote the density of those particles which are in single-occupancy sites at time t. The 
density of particles which are in pairs if then simply (p — psingie)- Note that the initial 
conditions selected earlier correspond to Psingie 

(0) = i. The single -occupancy property 
of a site is indicated by the value of £{2 — £) = 0, 1, 0, for £ = 0, 1, 2, respectively. This 
expression can also be written as a linear form in g and 9* (for k = 3). As a result 
we obtain the relation 



Psingie(t) = ^ {l - Re [(1 + iV3){aj+i{t)a;_,{t))] } (k = 3) (5) 

Let us point out that all the "linearity" properties also apply for k = 2. In 
fact, all calculations are then with real variables: the elements of Z2 are simply ±1. 
The dynamics (2) for the case k = 2 describes particles hopping on a dual lattice 
and annihilating pairwise on encounters, i.e., the cellular-automaton (simultaneous- 
dynamics) variant of the reaction A + A ^ inert; see [11]. Relation (4) is replaced 

by 



pit)=^[l-{aJ+,it)aJ-^m (k = 2) (6) 

Turning back to A; = 3, we note that the dynamical rules (2) imply 



{am{t)) = s = (7) 

for translationally invariant initial conditions. Since s = 0, we can work with the 
2-spin correlation function directly (for s one has to consider the "connected" 
part; the expressions are not much more complicated though). Thus, we define 



Gn{t) = {<Tj{t)a*+,{t)) (8) 

which does not depend on j due to translational invariance. This correlation function 
can be considered for n = 0, 2, 4, . . . only. Indeed, we have G_„(t) = G* (t) generally. 

By using (2), one can then derive the discrete-space and time diffusion equation 

for Gn{t), 



Gr^it + 1) - Gn{t) = - [Gn+2{t) - 2G,(t) + Gn-2(t)] (9) 
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which must be solved for n = 2, 4, . . . and t — 1,2,..., subject to the initial and 
boundary conditions 



Go{t>0) = l and Gn>o{t = 0) = (10) 

In writing (10), we used s = which implies a certain simplification of the prob- 
lem. Indeed, by (9) the correlation function Gn{t) is then real for these initial condi- 
tions. The form of the solution is thus essentially identical to that for k = 2. Relations 
(4), (5), (6) yield, for real k = 3 correlation function, 

[Pit)]k=3 = 2 [pW]fe=2 = 1 - G^it) (s = 0) (11) 



Psingle(t) = ^p(t) (A; = 3,s = 0) (12) 

The discrete diffusion equation (9) can be analyzed be generating function tech- 
niques [11]. The result for the density {k = 3) is 

Pit) = ^^i±^ (13) 
22^+1 [(t + l)!]2 ^ ^ 

This function decreases smoothly from the initial value 1, and in fact provides analytic 
continuation to all positive times t. The large-time asymptotic form is 



Pit » 1) ^ ^ (14) 



The fact that the density for A; = 3 turns out nearly identical to that for k = 
2 may at first seem disappointing. However, one must recall that all the various 
exact calculations for model Id reactions, based on the solution methods which involve 
evaluation of 2-spin correlation functions, always end up with continuous or discrete 



versions of the diffusion equation, — property shared with the Glauber Ising model in 
Id. Thus, the diversity of the functional forms obtained is limited [16]. The reaction 
scheme for = 3 is quite different from that ior k = 2. Indeed, if we denote particles 
in singly-occupied sites by A, and the pairs of particles, hopping as one unit, by B, 
then the reactions in the /c = 3 system are A + B ^ inert, A + A^B,B + B^A. 
Generally, such reaction scheme differs from the k = 2 reaction, A + A ^ inert, 
both in its mean-field rate equation form and in the time dependence of densities and 
correlations. 

Our results essentially suggest that for symmetric initial conditions (i.e., equal 
concentrations | of species A and B), the densities of the k = 3 and k = 2 reaction 
schemes are proportional, and both are fluctuation-dominated, ~ 1/ y/i, as opposed to 
the rate-equation prediction, ~ for large times. 
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